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ABSTRACT 

We present interpretation of the oscillation spectrum of the early B-type star y Pegasi, in which 
both low order p/g and high-order g-modes are observed. Using amplitudes and phases of the 
photometric and radial velocity variations, we identify /constrain the mode degree, I, for all 
14 detected frequencies. Seismic models fitting two pulsational frequencies corresponding to 
the modes I — 0, pi and I -\,%\ were constructed. This set of models contains those which 
reproduce also the empirical values of the complex nonadiabatic parameter / associated to 
these two mode frequencies. Unfortunately, there are no models reproducing the values of / 
for both frequencies simultaneously, regardless of model atmospheres, opacity data, chemical 
mixture as well as opacity enhancement in the Z- and Deep Opacity Bumps. Most probably, 
some modifications of the opacities in stellar interiors are still required. 
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1 INTRODUCTION 

The simultaneous presence of both acoustic (p) and internal grav- 
ity (g) modes in a star, allows potentially for probing almost 
the entire interior. Such hybrid oscillations have been detected 
in a few main sequence stars of early B spectral type: v Eridan i 
l lHandler et alj|2004 lAerts et alj|2004 Ijerzvkiewicz et al.ll2005h . 



12 Lacertae ^Handler et all I2006J: iDesmet et alJ l2009h . y Pegasi 
(Handler e t"aij|2009h . HD50230 (Degroote eTalJl2012h . HD43317 
(Papics et al. 2012) and many candidates have been found (Pigul- 
ski& Pojmanski 2008, Degroote et al. 2009, Balona et al. 2011). It 
was only possible thanks to multisite campai gns and space obser - 
vations from the satellite m issions like MOST (Walker et al. 2003), 
CoRoT teaglin et alj|2006l) and Kepler <Koch et alj|20ld) . The ex- 
istence of both low order p/g modes as well as high-order g-modes 
makes the early B-type pulsators the attractive targets for astero- 
seismic studies. 

The basic goal of asteroseismology is to find pulsational 
models with eigenfrequency reproducing the observational 
values. A more advanced seismic model should additionally 
account for instability and properties of oscillati on modes . 
Daszyn ska-Daszkiewicz. Dziembowski & Pamvatnvkri (2003J, 
20051) introduced a new asteroseismic tool associated with 
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each pulsational frequency. It is a ratio of the amplitude of the 
bolometric flux perturbations to the radial displacement, called 
the nonadiabatic /-parameter. The theoretical values of / are 
obtained from stellar pulsation computations and their empiri- 
cal counterparts are derived from multicolour photometry and 
radial velocity data simultaneously with the mode degree, I. 
The parallel fitting of pulsational frequencies and corresponding 
values of the /-parameter was termed comp lex asteroseimology by 
iDaszvnska-Daszkiewicz & Walczakl y009). The hybrid pulsators 
are of particular interest for such an in-depth modelling because 
the /-parameter depends on pulsational frequency and the shape 
of eigenfunctions. Therefore its behavior is very different for the 
different frequency ranges. In the case of low order p/g modes (the 
high frequency range), the values of/ are independent of the mode 
degree I and slowly vary with the frequency. On the contrary, for 
high order g-modes (the low frequency range), the values of / are 
strongly ^-dependent and change rapidly with the frequency. 

Complex seismic modeling has been already applied to 
the B Cep star 9 Oph (Daszynska-Daszkiewicz & Walczakl 
120091) and to the two /? C ep/SP B stars: v Eri 
I Daszvhska-Daszkiewicz & Walczakl l2010f) and 12 Lac 
(Daszynska-Daszkiewicz, Szewczuk & Walczak 2013). In 
the case of 6 Oph, we got a strong preference for the OPAL 
dlglesias & Rogers! 1 1996T) opacity tables. From the analysis of v 
Eri a contradictory result was obtained: the {S Cep-type modes 
indicated the OPAL da ta whereas the SPB-type modes preferred 
the OP dSeatonl 2005) opacities. For 12 Lac, somewhat better 
agreement was found with the OP data. It can indicate, that 
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the used opacity tables are not properly defined in the whole 
range of temperature and various physical conditions. Up to now, 
the sample of analyzed stars is too small to draw more general 
conclusions. 

In this paper, we present mode identification and complex 
seismic modeling of 7 Pegasi. The an alysis of the space based 
observations from the MOST satellite JHandler et"aj] |2009|) and 
ground based photometry and spectroscopy (Handler 2009) led to 
the discovery of 14 pulsational frequencies with 8 typical for the 
p Cephei stars and 6 typical for the Slowly Pulsating B-type stars 
(SPB). Because the star is a very slow rotator, here, we neglect all 
effects of rotation on pulsations. 

In Section 2, we give a short description of the star. Section 3 
contains results on mode identification for all detected pulsational 
frequencies using two approaches. Results of our seismic modeling 
are presented in Section 4. Conclusions are summarised in Section 
5. 



2 THE HYBRID PULSATOR: 7 PEGASI 

7 Peg (HR 39, HD 886) is a pulsating star of B2IV spectral type 
with the brightness of V=2.83 mag. The radial velocity varia- 
tions were discovered bv lBurnsl ( ll91ll) but a very small amplitude 
made the variabil ity uncertain. Series of spectrograms obtained by 
iMcNamaral ( 119530 showed, that this star is a variable with the pe- 
riod of about 3 h 38 m and the am p litude of the radial velocity vari- 
ations of 3.5 km/s. IMcNamaral dl953l) attributed 7 Peg to the yS 
Cephei class of variable stars on the basis of its short period and 
spectral type. The follow-up spectroscopic and photometric obser- 
vations confirmed changes of the rad ial velocity (IMcNamara 1955, 
1 19561 : ISandberg & McNamaral Il960h and led to the discovery of 
the li ght variations i n the yellow filter w ith the range of 0.015 
mag (Williams 1 9541 ). I Jerzvkiewica d 19701) determined the follow- 
ing light ranges in the UB V passband s: 0.027±0.001, 0.018±0.001 
and 0. 017 ± 0.002 mag, respectively. ISarevan. Valtier & Le Contell 
d 1975b determined period more precisely, P = 0.1 517501(3) d, and 
confirm ed i ts constancy within 0.06 s per century. ISmith & McCalll 
dl978l) and ICugier. Dziembowski & Pamyatnykh! d 19941) showed, 
that the variability of 7 Peg can be explained by the radial pulsation. 
In the catalogue of the galactic /3 Cephei stars by 



IStankov & Handled ( I2005I) , 7 Peg is the lowest mass variable, close 
to the Slowly Pulsating B stars (SPB) instability strip. 7 Peg had 
been thought to be o ne of the monoperiodic /3 Cep star until 
IChapellier et al.l (2006) reported three additional pulsational fre- 
quencie s; one in the 13 Cephei r ange and two in the SPB frequency 
domain. IChapellier et al.l ( 120061) claimed also that the star is a spec- 
troscopic binary with an orbital period of 370.5 d. However, the 
analysis of data collected by the MOST satellite and from ground- 
based photometric^ and spectroscopic observations ( Handler et al.l 
2009; Handler 2009), showed that 7 Peg is a single star and the hy- 
pothetical orbital variations can be explained by the high-order g- 
mode pulsat ion. Moreo ver, Handler et al. (2009) confirmed the fre- 
quencies of IChapellier et al . (2006) and discovered ten new ones. 
Consequently, we know that the light variations of 7 Peg are caused 
by at least 14 pulsational frequencies: 8 of them are of the /? Cep 
type and 6 are of the SPB type. Thus, 7 Peg is the hybrid pulsator 
in which low-order pressure/mixed modes and high-order gravity 
modes are excited simultaneously. 

The ro t ational velocity of 7 Peg was determined by 
iTelting et al.l d2006l) who derived V rot sin i ts from the Si- 
III lines. Similar result was given by Pandey et al. (2011) and 
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Figure 1, The observational error boxes of y Peg in the HR diagram. The 
evolutionary tracks were computed for metallicity of Z = 0.015 and two 
values of the overshooting parameter, a m . Lines of the constant period 
(0.15175 d) for the fundamental (n = 1) and first overtone (n = 2) radial 
mode are also drawn. Points labeled as a diamond and asterisk are discussed 
in Section 4. 



iNieva & Przvbillal d2012l). who got V rot = 6 km/s and V rol = 9 ± 2 
km/s, respectively. Handler et alj J2009I) assuming that two pulsa- 
tional frequencies are components of the t = 1, gi triplet, deter- 
mined V rot ~ 3 km/s from the rotational splitting. 

Det erminations of the heavy elements abundance are am- 
biguous. iDaszvhska-Daszkiewicz & Niemczura I (120051) obtained 
the metal abundance of [m/H]= -0.04 ± 0.08 (equivalent to the 
metallicity parameter ofZ « 0.017 ± 0.004) from the IUE ultra- 
violet spectra. Fr om the optical spectra, iMorel et al.l d200 6) and 
IPandev et al] d201ll) derived Z = 0.009 + 0.002 and Z = 0.010 ± 
0.002, respectively. The recent determinations fro m the optical and 
ultraviolet IUE spectra give Z = 0.017 + 0-001 dWu et al.|[2oTTb 
and Z = 0.014 ± 0.002 dNieva & Przvbillall2012ri , while the anal- 

i sis of the ultraviolet HS T spectra indicates Z = 0.016 ± 0.003 
Koleva & Vazdekisll2012l) . 

In FigJj] we show the observational error box of y Peg in 
the HR diagram. We i ncluded the mos t recent det ermination of ef- 
fect ive temperature by ;Wu et al. ( 201ll) , lKoleva & Vazdekisld2012r) 
and lNieva & Przvbillal d2012D . The luminosit y was calculated from 
the Hipparcos parallax, n = 8.33 ± 0.53 mas dvan Leeuwenll2007l) . 
and the bolometric correction from Flower (1996). The total er- 
ror box is as follows: log T c9 = 4.325 ± 0.026 and log L/L = 
3.744 + 0.090. 

We depicted also the evolutionary tracks from the zero- 
age main sequence (ZAMS) to the terminal-age main sequence 
(TAMS) for masses M = 8.0, 8.5, 9.0 and 9.5M Q , the hydro- 
gen abundance of X = 0.7, metallicity of Z = 0.015, two val- 
ues of the overshooting parameter, ff ov = 0.0 and 0.3, the initial 
equatorial rotational velocity of V rot = 3 km/s and the element 
mixture by Asplund et al. (2009), hereafter AGSS09. The tracks 
were computed by means of the Warsaw-New Jersey evolutionary 
code ( e.g Pamyatnyk h et alj|l998f) adopting the OPAL opacities 
dlglesias & Rogers! 1 19961) . Lines labeled as n = 1 and n = 2 and 
models marked with a diamond and asterisk will be discussed later 
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Table 1. Pulsational observables of y Peg. In the first column, we give the pulsational frequencies of y Peg. In the following columns there are values of the 
Stromgren u, v, y amplitudes and phases. In the last two columns there is the radial velocity amplitude and corresponding phase. 



frequency 


A u 


<Pu 


A v 


<Pv 


Ay 


tPy 


Avrad 


fVrad 


[c/d] 


[mmag] 


[rad] 


[mmag] 


[rad] 


[mmag] 


[rad] 


[km/s] 


[rad] 


v, =6.58974(2) 


12.26(13) 


2.669(11) 


6.74(10) 


2.628(15) 


6.074(94) 


2.607(15) 


3.3582(69) 


4.490(2) 


v 2 =0.63551(10) 


2.44(13) 


0.870(57) 


1.55(10) 


0.734(67) 


1.269(92) 


0.947(76) 


0.4992(70) 


4.735(14) 


v 3 =0.68241(7) 


2.21(13) 


5.607(61) 


1.64(10) 


5.681(63) 


1.500(94) 


5.644(63) 


0.7364(72) 


3.402(10) 


v 4 =0.73940(10) 


1.77(13) 


3.542(75) 


1.18(10) 


3.604(89) 


1.166(95) 


3.782(83) 


0.5214(71) 


1.468(14) 


vj=6.01616(14) 


1.58(14) 


4.592(87) 


1.12(10) 


4.664(92) 


0.828(97) 


4.54(12) 


0.3612(85) 


0.148(24) 


v 6 =0.88550(7) 


1.54(14) 


3.733(87) 


1.10(10) 


4.014(94) 


0.896(94) 


4.00(11) 


0.7236(68) 


0.873(10) 


v 7 =6.9776(5) 


0.48(13) 


3.08(28) 


0.25(10) 


3.27(42) 


0.328(95) 


3.36(29) 


0.0962(72) 


4.930(75) 


v 8 =0.91442(ll) 


0.94(13) 


4.45(15) 


0.60(10) 


4.80(17) 


0.417(95) 


4.95(23) 


0.4652(68) 


1.191(16) 


vq =6.5150(8) 


0.18(13) 


3.46(73) 


0.15(10) 


2.44(68) 


0.258(94) 


2.64(37) 


0.0628(70) 


4.15(11) 


v 10 =8.1861(8) 


0.26(13) 


4.92(50) 


0.26(10) 


4.89(40) 


0.143(94) 


5.57(66) 


0.0639(69) 


1.60(11) 


vi i =0.8352(3) 


0.82(13) 


1.30(16) 


0.60(10) 


0.97(17) 


0.700(93) 


1.43(14) 


0.1785(67) 


4.743(41) 


vi2=6.0273(5) 


0.41(14) 


1.55(34) 


0.35(10) 


0.53(30) 


0.389(95) 


0.57(25) 


0.1104(82) 


1.823(75) 


Vi3=9. 1092(12) 


0.38(13) 


3.61(34) 


0.26(10) 


3.38(39) 


0.279(94) 


3.24(34) 


0.0408(68) 


4.51(17) 


vi4=8.552(2) 


0.20(13) 


0.06(67) 


0.37(10) 


0.59(28) 


0.222(94) 


5.54(42) 


0.0260(69) 


1.70(27) 



3 IDENTIFICATION OF OSCILLATION MODES 

To identify the mode degree, I, we made use of the light variations 
in the Stromgre n uvy passbands I Handlenl2009B and the radial ve- 
locity changes ( Handler et alj2009l) . In TableQ] we give amplitudes 
and phases of the light and radial velocity variations determined by 
fitting 14 frequencies detected in the MOST data. Therefore, our 
values of the amplitudes and phases can differe from those deter- 
mined by Handler (2009). 

We performed identification of the degree, 6, by applying 
two methods. In the first case, we compared theoretical and 
observational values of the amplitude ratios and phase differ- 
ences between the available passbands and relied on the theo- 
retical values of the /-parameter. In the second approach, we 
used amplitudes and phases themselves and the /-parameter 
was deter mined from the observations together with the mode 
degree, I i Daszvnska-Daszkiewicz. Dziembowski & Pamvatnvkhl 



2003, 2005). In the case of B-type pulsators, the latter method de 
mands the radial velocity measurements to get a unique identifica- 
tion of I. 

In both cases we need input f rom mod el atmospheres. Here, 
we present results obtained with the lKuruczl ( |2004[) models and the 
microturbulent velocity of £, = 2 km/s. Coefficients of the non- 
linear limb darkening law were adopted from Claret (2000). Iden- 
tification of (, does not change if other val ue of f ( or the non-LTE 
model atmospheres dLanz & Hubeny|2007r) were used. The theoret- 
ical values of the / -parameter are calcul ated with the nonadiabatic 
pulsational code of lDziembowskil J19771) . 

In Table[2] we give the most probable values of the mode de- 
grees for the y Peg frequencies from the two approaches. The dom- 
inant mode, Vi, is certainly radial. Also the identification of the 
V3, vs and V(, is unique; they are modes with £=1,1 and 2, respec- 
tively. The frequencies v 2 , v 4 and V\ \ can be either dipoles (C = 1) or 
quadruples (C = 2). We got also two possibilities for the frequencies 
vg and Vi2. They can be £ = 3 or 2 and £ = 2 or 5 modes, respec- 
tively. Identification of the remaining frequencies are ambiguous. 
In the case of V7, vg, v^ and V13 we can say only that their mode 
degrees should be different from 4 and 6, and in the case of Vu 
only degrees higher than 3 are possible. It is interesting that Vu is 
close to the first overtone radial mode in many 7 Peg models, but 
the degree I = was excluded by our photometric identification. 



Table 2. Identification of £ from the two methods. In the first column, we 
give the pulsational frequencies of y Peg, the second and third columns con- 
tain the mode degrees, t, identified from photometric observables using the 
theoretical values of / and from photometric observables and radial veloc- 
ity data using the empirical values of /, respectively. In the last column, we 
put values of t consistent for both approaches. 



frequency 


phot. 


phot.+V r ad 


e 


[c/d] 


theoret. / 


empir. / 


v, =6.58974 


£=0 


£=0 





v 2 =0.63551 


£=1,2,4 


£=2,1,3,5 


2,1 


v 3 =0.68241 


e=\ 


£=1 


1 


v 4 =0.73940 


£=1,2 


£=1,2 


1.2 


v 5 =6.01616 


C=\ 


£=1,3 


1 


v 6 =0.88550 


1=2 


£=2 


2 


v 7 =6.9776 


1=1 


£*4,6 


*4,6 


y 8 =0.91442 


£^4 


£=3,2,5 


3,2 


v 9 =6.5150 


1=1 


£*4,6 


*4,6 


vio=8.1861 


1=1 


£*4,6 


*4,6 


vi 1 =0.8352 


£=1,2 


£=1,2 


1,2 


vi2=6.0273 


£=? 


£=2,5 


2,5 


vi3=9.1092 


1=1 


£*4,6 


*4,6 


v 14 =8.552 


£>4 


£=? 


5*4 



In the HR diagram (Fig.Q~|l, we plot lines of a constant pe- 
riod (0.15175 d) corresponding to v\ for the fundamental and first 
overtone radial modes. As we can see the pi line is included only 
marginally in the error box. Nevertheless, two possibilities have 
to be considered. The discrimination of the radial order, n, can be 
done by a comparison of the empirical and theoretical values of 
the /-parameter (Daszyhska-Daszkiewicz, Dziembowski & Pamy- 
atnykh 2005). In Fig.[2] we plotted this comparison on the complex 
plane assuming that the dominant frequency is the fundamental (the 
top panel) or the first overtone mode (the bottom panel). The theo- 
retical values of / were calculated for three different metallicities, 
Z = 0.010, 0.015, and 0.020, overshooting parameter a m = 0.0, 
hydrogen abundance of X = 0.7 and the OPAL data. In all cases, 
we included models that are inside the observational error box of 7 
Peg. In the upper panel, we showed the line of the constant insta- 
billity parameter 77, defined as (Stellingwerf 1978): 
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Figure 2. Comparison of empirical and theoretical values of / for the dom- 
inant frequency, v\, on the complex plane (/r, fi) for models from inside 
the total error box. The OPAL data and AGSS09 mixture were adopted. The 
top and bottom panels correspond to the hypothesis of the fundamental and 
first overtone radial mode, respectively. In the upper panel, models to the 
right of the r\ = line are stable; in the lower panel all models are stable. 
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Figure 3. The overshooting parameter, a ov , as a function of metallicity, Z, 
for seismic models of y Peg found from the fitting of the V[ frequency (the 
t = 0, pi mode) and the v^ frequency (the C = 1, gi mode) for the hydrogen 
abundance X = 0.7 and OPAL opacities. We plot also the lines of constant 
mass, effective temperature and luminosity. Indicated vales of log r c (f and 
log L/Lq are for the center and 1 <x error of the total error box. Models below 
the r\ = lines are unstable. 



tended by a requirement of reproducing also the empirical values 
of the nonadiabatic parameter /, which corresponds to each mode 
frequency. 

In all comput ations, we used the Wars aw-New Jersey stellar 
evolution code (e.g. Pamyatnykh et al. 1998) and nonadiabatic pul- 
sational code of iDziem bowski ( 1977). Seismic models were cal- 
culated with both the OPAL and OP data. We assumed the initial 
rotational velocity of 3 km/s a nd two chemical mixture s: AGSS09 
and that determined for y Peg ( Nieva & Przvbillall2012i) . hereafter 
the y Peg mixture. Due to the excitation problem of some modes, 
we tested also the opacity enhancement near the Z-bump and Deep 
Opacity Bump (DOB). 



IV 



r 



dW 

th- 



eir 



(1) 



where W is the work integral and R is the stellar radius. Models lo- 
cated to the right of the 77 = line are unstable, whereas in the lower 
panel all models are stable. As we can see, agreement between the 
empirical and theoretical values of the /-parameter can be achieved 
only if Vi is the radial fundamental mode, n = 1. Moreover, con- 
straints on the metal abundance were obtained, viz Z e (0.012, 
0.014). The conclusion about the radial order of V\ did not change 
if the OP opacities were used. In this case, we got the metallicity 
in the range of Z € (0.010, 0.012). Moreover, identification of the 
radial order is independent of the value of the core overshooting pa- 
rameter, a m . For the higher values of a m , the allowed metallicity 
range is shifted to the smaller value of Z. 



4 COMPLEX ASTEROSEISMOLOGY 

Basic seismic modeling consists in fitting pulsational frequency 
taking into account instability condition. This approach can be ex- 



4.1 Fitting the centroid frequencies of y Peg 

In the first step we searched for models fitting the value of the dom- 
inant frequency, Vi, identified as the radial fundamental mode in 
previous section. Next, from this set of models we selected only 
those, which fit also the observational value of the frequency V5. 
In the whole set of models, we found that if V\ is the radial funda- 
mental mode, then V5 can be only the dipole gi mode. This identifi- 
cation of V5 is unique, because the difference between frequencies 
of consecutive dipole modes is of the order of 1 c/d. Identification 
of the radial order is also independent of the azimuthal number, m, 
since the rotational splitting is of the order of 0.005 c/d. The fre- 
quencies V\ and v 5 were chosen because they are well identified 
low-order pressure/mixed modes. High-order gravity modes have a 
very dense frequency spectrum and fixing the radial order is prob- 
lematic. The accuracy of the fitting V\ and v 5 is of the order of 10~ 5 
c/d, which is equivalent to the observational errors. 

These seismic models of y Peg, calculated with the hydro- 
gen abundance of X = 0.7, the AGSS09 chemical composition and 
OPAL opacities, are shown in Fig.[3]on the a m vs. Z plane. We 
depicted the lines of constant mass, effective temperature and lu- 
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minosity, as well as the instability borders for the radial mode, V! 
(thick solid line) and for the dipole mode, v 5 (thick dashed line). We 
indicated masses from 6.8 up to 9.8M Q and the values of logT e (j- 
and log L/L corresponding to the center and edges of the total 
error box. Models that are located below the 77 = lines are un- 
stable. Almost all models between the lines log L/L Q = 3.654 and 
log L/L Q = 3.833 are inside the observational error box of 7 Peg. 

As we can see, for a given mass, models with higher metal- 
licity require smaller values of the overshooting parameter, effec- 
tive temperature and luminosity. The same results were obtained 
for 6 Oph (Daszyriska-Daszkiewicz & Walczak 2009) and 12 Lac 
(Daszyhska-Daszkiewicz, Szewczuk & Walczak 2013). Lines of 
constant effective temperature and luminosity are nearly parallel to 
each other and more steep than lines of constant mass. Having such 
a huge number of seismic models fitting two frequencies of 7 Peg, 
we determined an approximate relation between their parameters: 



a m = -21.36(33)Z - 0.2613(38)M + 2.623(30). 



(2) 



The very important point of the seismic analysis is mode in- 
stability. As can be seen in Fig.[3] all models with Z < 0.01 have 
the modes £ = 0, pi and £ = 1, gi stable. It is caused by a small 
abundance of heavy elements, which are crucial for exciting pulsa- 
tions in the B-type stars. A lot of unstable models appear for higher 
value of metallicity, Z > 0.01, and overshooting less than about 
0.35 for the radial mode vi and about 0.48 for the dipole mode V5. 
Models with high overshooting are stable because they have too 
low masses and fall outside the instability region of the /3 Cep stars 
( IPamvatnvkhll 1 9991) . For more thorough dis cussion of the mode in- 
stability analysis see lZdravkov et al.l ( 120131) . 

In Fig.|4]and|5] we plotted the same as in Fig.[3] but we marked 
also lines of models fitting some of seven frequencies: v 2 , v 3 , v 4 , 
v 6, v 8> y n and V13. In most cases, we assumed that these modes are 
axisymmetric (m = 0) because there is no determination of m. For 
v 3 we considered two possibilities, m = an d m = 1, since ther e 
are some indications that the mode is prograde (IPandev et al.l201 ll) . 
We considered simple rotational splitting in order to calculate the 
centroid frequency, that is v„ e „, = v n( + m{\ - C„f)v Iot , where C„t is 
the Ledoux constant and v rot is the rotational frequency. 

Models located on each line fit simultaneously three frequen- 
cies. Fig.|4] includes frequencies which were identified as £ = 1 
or 3, whereas Fig.|5]- those which have the degree I = 2. These 
seven frequencies were fitted with an accuracy of 10~ 4 - 10~ 3 c/d, 
equivalent to the observational errors. 

In the case of high-order g-modes (vt, V3, v 4 , vg, vg, vn), we 
usually had to include more than one radial order. In a considered 
range of metallicity and overshooting parameters, the frequency V3 
is the I = 1, gn mode (in both cases, for m = and m = 1) and 
the frequency v 6 is the £ = 2, gi6 or g 17 mode. Identification of 
the mode degree of vi, v 4 , vg and Vn is ambiguous and we had to 
consider more possibilities. For V2 we got the modes £ = 1, gi2 and 
1 = 2, g 2 2, g23> g24> f° r y u - the modes I = 1, g 9 and 1 = 2, g 17 - 
The frequency vg can be a mode with £ = 2 or (. = 3. In the former 
case, it has to be the gi5 or gi6 mode and in the later - the g2i or 
g 2 2 mode. If the frequency v 4 is a dipole, then it has to be gi , but 
if it is a quadruple, it is g^ or g2o- Because in many models the 
pressure mode C = 1, P2 has a frequency close to V13, we plotted 
also a model line reproducing this frequency. 

Models fitting three frequencies require higher overshooting, 
effective temperature, luminosity and mass if metallicity decrease. 
Unfortunately, high-order g-modes are usually stable. The modes 
V2{t = 2, g23) and V(,{t = 2, gi6) are exited if the metallicity param- 



eter Z > 0.02; v 2 (l = 2, g 24 ), v 4 (£ = 2, g 20 ) and Vi (£ = 3, g 2 i) if 
Z > 0.016 and v 6 (£ = 2, g 17 ) and v s (£ = 3, g 22 ) if Z > 0.013. 

Moreover, we were able to find a lot of models fitting four 
frequencies and two models fitting five frequencies. A model with 
M = 8.093M Q , log r cff = 4.3286, a m = 0.27 and Z = 0.01 1 fits the 
modes v t (£ = 0,pi), v 5 (i = 1, gO, v 3 (I = 1, g n , m = 0), v 6 {I = 2, 
gi6) and vn (( = 2, g i7 ). We will call it MODEL1 and it is marked 
with a diamond in Fig.fj] The parameters of the other model fitting 
five modes are M = 7.895M , logr cff = 4.31956, a m = 0.30 and 
Z = 0.012. The fitted modes are: v\{€ = 0, pi), v 5 (£ = 1, gi), v 2 
(£ = 1, gi2), vn {£ = 1, gg) and V13 (C = 1, P2). It is worth to note 
that these two models are located inside the observational error box 
of 7 Peg. However, neither of these two models can reproduce the 
frequency v 4 , what can indicate that it is not an axisymmetric mode, 
as we have assumed. 

In Fig.[6] we present the detailed comparison of the observa- 
tional frequency spectrum and their possible theoretical counter- 
parts for MODEL1 described above. We considered modes with C 
from to 3. Each theoretical frequency peak is labeled with the 
radial order, n. The well identified frequencies have the best theo- 
retical counterparts for mode degrees generally consistent with the 
photometric identification. The exception is V2. Despite the fact that 
it was identified as a dipole or quadruple mode, the best match is 
for the ( = 3 mode. A comparison with the theoretical spectrum 
may suggest, that this is the prograde mode (m > 0). The frequency 
V12 can be the prograde mode, I = \, m=-l, gi as assumed by Han- 
dler et al. (2009). The frequencies V9 and v^ seems to be the 1 = 2, 
gi and po modes, respectively. The frequencies v 7 and V13 can be 
consecutive dipole modes, i.e., pi and P2. We can see also, that the 
frequency v 14 is consistent with the first overtone radial mode, but, 
as we mentioned in Section[3] it does not agree with our identifica- 
tion. The other close mode to v [4 is I = 6, gi (not shown in Fig. 6), 
which is not inconsistent with photometric identification. 

Let us now check the instability conditions for MODEL 1. In 
the left panel of Fig.[7] we plotted the instability parameter, rj, as 
a function of frequency for modes with £ = — 4. For a compar- 
ison, in the right panel of Fig.|7] we plotted the same for a model 
with very similar parameters but computed with the OP data. Short, 
vertical lines represent the observed frequencies of 7 Peg. Modes 
with 77 > are excited in the model. As we can see, with the OPAL 
opacities only frequencies in the range from about 5 up to 7 c/d are 
unstable. In the domain of the high order g-modes only high-degree 
modes (£ ^ 4) are excited with frequencies v > 1 c/d. The frequen- 
cies in the range 8-9.5 c/d and 0.6-0.9 c/d are stable. Slightly better 
situation appears with the OP data. The model computed with these 
opacity data excite high-order g-modes with £ > 2. But neither with 
the OP nor OPAL date we were able to excite modes with frequen- 
cies higher than about 7 c/d. It is tempting to explain this unsatis- 
factory results assuming higher opacities in the driving zones, but, 
as we will discuss it in the Section |4~3l it causes other problems. 



4.2 Fitting the nonadiabatic parameter, / 

As we have mentioned in the Introduction, an adequate seismic 
model should reproduce not only pulsational frequencies but also 
other characteristic parameters. Such parameter, which is deter- 
minable from observations, is the /-parameter associated with each 
frequency. The theoretical values of / depend mostly on properties 
of the driving zone. 

To derive the empirical values of the /-parameter, we used the 
LTE models of stellar atmospheres with the microturbulent velocity 
of i;t — 2 km/s. The results do not differ qualitatively from those 
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Figure 4. The same as in Fig.[5] but we added lines of seismic models fitting the high-order g-mode frequencies which were identified as t = 1 or ( = 3 modes. 
These are: V2 (A), V3 (<i), V4 (V), vg (♦) and vn (□). The line of models fitting the p-mode frequency V13 (f = 1, P2) is plotted as (V). The grey area indicates 
models located inside the total error box of y Peg. In the case of V3 two values of m were considered: m = and m = 1. 



comp uted with the non-LTE atmosphere models dLanz & Hubenvl 
2007) as well as with higher microturbulent velocity £j, = 8 km/s. 

In Fig.[8] we put seismic models fitting, within the observa- 
tional errors, the empirical values of the /-parameter (hatched ar- 
eas) for the radial fundamental mode (labeled as /(vi )) and for the 
dipole gi mode (labeled as /(V5)). The area /(V5) is larger than 
/(vi) because of larger observational errors in the photometric am- 
plitudes and phases of the frequency V5. The left and right panels 
of Fig.[8] correspond to computations with the OPAL and OP ta- 
bles, respectively. As we can see, the requirement of fitting the /- 
parameter for the v, mode reduce significantly the allowed range of 
stellar parameters of 7 Peg. With the OPAL tables and the hydro- 
gen abundance X = 0.7, only models with a m « 0.3 and Z as 0.013 
can be considered. The OPAL models fitting the /-parameter for 
the V[ mode are inside the observational error box of 7 Peg, while 
all models fitting the /-parameter for the V5 mode are outside the 
box. In the case of the OP models, the areas /(vi) and part of f(v 5 ) 
are inside the error box. 

Results for the OPAL and OP tables are similar and the biggest 
difference is in the position of the instability borders. For the OP 
models, these lines are shifted to the lower values of metallicity, Z, 
and the overshooting parameter, a m . Furthermore, for a given value 
of Z and a m , models with the OP tables have slightly higher mass 
whereas effective temperature and luminosity are almost the same. 
Moreover, the areas indicating models fitting the /-parameter ap- 
pear at the lower overshooting and metallicity. Unfortunately, there 



is no model fitting the /-parameter for both v t and v 5 , simultane- 
ously; neither with the OPAL nor OP data. 

In the case of the most high-order g-modes, the empirical 
and theoretical values of the real part of / agreed and the imag- 
inary part differed significantly. This occurs both with the OPAL 
and OP opacities. Only for a mode V(, we were able to find some 
seismic models which fit its empirical value of / (both with the 
OPAL and OP data). These models, obtained with the OP data, have 
a m » 0.31 and Z ss 0.008 - 0.009, and are located inside the obser- 
vational error box of 7 Peg. However, they do not overlap neither 
with the area f(v\) nor /(vj). The OPAL models reproducing the 
/-parameter of V(, are outside the allowed range of parameters and 
have a m ss 0.51 and Z « 0.010 - 0.012. 

In the next step, we checked the effects of the hydrogen 
abundance and the heavy elements composition. In comparison 
with the results computed with the standard hydrogen abundance, 
X = 0.7, models with X = 0.75 for a given value of Z and a m 
have larger masses (about 0.3M o ), smaller effective temperatures 
(AlogT e ff ~ 0.18) and smaller luminosities (AlogL/L G ~ 0.04). 
Also models fitting the /-parameters required much smaller value 
of the overshooting parameter. Still, there were no models fitting 
the /-parameter for both V\ and V5 frequency simultaneously. 

Then, we changed the chemical mixture. In Fig.[9] we show 
seismic models c alculated with the chemic al composition of 7 Peg 
as determined bv lNieva & Przvbillal ( 12012b . The 7 Peg mixture has 
smaller abundance of iron-group elements than the Sun. These el- 
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Figure 5. The same as in Fig.[4] but models fitting frequencies identified as high-order g-modes with C = 2 were added. These are: V2 (A), v\ (►), vg (ir), v% 
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Figure 6. Comparison of the observed frequencies of y Peg (upper panel) with the theoretical counterparts (bottom panel) corresponding to MODEL 1 of y 
Peg which has the following parameters: M = 8.093M Q , log T c g = 4.3286, a m = 0.27 and Z = 0.01 1. The height of lines in the upper panel corresponds to 
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Figure 7. Run of the instability parameter, r\, with frequency for MODEL 1 (the left panel) and a model with very similar parameters but computed with the 
OP tables (the right panel). The vertical lines corresponds to the observed oscillation spectrum of y Peg. 
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Figure 8. The same as in Fig.[5]but we marked models fitting the empirical values of the nonadiabatic /-parameter (hatched areas) of the radial fundamental 
mode (labeled as f(y\)) and of the dipole mode g[ (labeled as /(V5)). In the left panel we used the OPAL data and in the right panel the OP data. 



ements are extremely important for exciting pulsation and because 
of their deficiency, the instability region is smaller. Nevertheless, 
there are still a lot of unstable models fitting two pulsational fre- 
quencies and the /-parameter corresponding to V\ or v 5 . In this case 
seismic models fitting the values of the /-parameter require slightly 
less efficient core overshooting. 

Because the formal errors of / can be underestimated, in 
Fig.|9] we marked additionally models fitting the /-parameters 
within 2cr, (shaded-in, dark areas around the hatched regions) and 
3cr (shaded-in, bright areas around the hatched regions), where a is 
the empirical error of /. Now, the majority of models reproducing 
/ for vi fit also, within 3<r, the /-parameter corresponding to V5. 

It is important to add, that in all considered cases, we were 
able to find also a model fitting, within lcr error, the /- parame- 
ter for Vi and the real part of / for V5. The disagreement remains 



for the imaginary part of /. This is clearly visible in Fig.[lO] where 
we compare the theoretical and empirical values of / in a func- 
tion of the frequency. The left and right panels show the real (/r) 
and imaginary (/j) part, respectively. The parameters of the cho- 
sen seismic model marked with an asterisk in Fig.[T] and [9] are: 
M = 7.968M Q , logr cff = 4.3134, logL/L Q = 3.6728, X = 0.71, 
metallicity Z = 0.0135, a m = 0.25, the OPAL data and 7 Peg mix- 
ture. We called it MODEL2. Modes with the degrees (, = 0, 1, 2, 3, 4 
were considered. 

We marked the empirical values of the /-parameters for V| , V5 
and six high-order g-modes: v 2 (£ = 2 and 1), v 3 (t = 1), 1/4 ({ = 2 
and 1), v 6 (l = 2), v 8 (t = 2 and 3), v u (I = 1 and 2). The value of 
the imaginary part of the empirical /-parameter for v 8 it - 3) is out 
of the scale. The empirical values of / corresponding to different 
degrees have different symbols. We can conclude that only in the 
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Figure 9. The same as in the left panel of Fig.[8]but for the chemical com- 
position of y Peg as determined by Nieva & Przvbilla (2012). The dark 
shaded-in and bright shaded-in areas around the hatched regions indicate 
models fitting the /-parameters within 2<r and 3<r, respectively. The aster- 
isk marks model chosen for a comparison of the values of /. 



case of the real part of / an agreement between the theoretical and 
em pirical values of / is quite good. The same result was obtained 
by Daszyhska-Daszkiewicz, Dziembowski & Pamyatnykh (2005) 
and lDaszvhska-Daszkiewicz & Walczakl 120101) for v Eri. 

The values of the /-parameter for p modes are almost indepen- 
dent of the mode degree, I, and change slowly with the mode fre- 
quency, while for high-order g-modes /-parameter depend strongly 
on t and changes rapidly with v. The fact, that the empirical values 
of the real part of /, especially for high-order g-modes, are located 
almost perfectly along the lines of theoretical counterparts are very 
encouraging for further studies and it shows a great potential of the 
/-parameter, particularly for the SPB-type modes. 

The empirical values of the nonadiabatic /-parameter, plotted 
in Fig.[lO] are listed in Table[3] In the fifth column we give also 
the values of the intrinsic amplitude, |s|, multiplied by spherical 
harmonic, Yf{i, 0), where i is the inclination angle. In the case of 
the radial mode, V\ , we have the exact value of |e|, whereas for other 
modes we can get only their lower limit. As we can see, the intrinsic 
amplitude of the dominant radial mode is very small and amounts to 
~0.26 per cent of the stellar radius and it is about three times larger 
than for other modes. The last column contains the values of the 
discriminant, %%, which measures goodness of the fit between the 
theoretical and observational values of the photometric amplitudes 
and phases. 



4.3 Effect of the opacity enhancement 

In the previous section, we have tried to find seismic models of y 
Peg reproducing two frequencies and corresponding values of the 
/-parameter. Using the standard opacity data and \a error in / we 
did not succeed in finding a model fitting the values of / for v\ and 
vj simultaneously. Because of this, we decided to change artificially 
the OPAL opacity tables. Two cases were examined. 

In the first case we tested the 50% opacity enhancement near 
the Z-bump, i.e., at log T ~ 5.3 - 5.5. The lines of constant ef- 



Table 3. The empirical values of the real and imaginary part of /, and the 
absolute value of the intrinsic mode amplitude |e| for 8 frequencies of y 
Peg. Columns from left to right are: pulsational frequency, mode degree, 
I, empirical values of the real, /r, and imaginary, f\, part of the nonadia- 
batic parameter, the intrinsic amplitude multiplied by spherical harmonic, 
|e| = \sY"'(i, 0)|, and the discriminator, ^ 2 , describing the goodness of the 
fit. These quantities were derived for the parameters of MODEL2, marked 
with asterisk in Fig.[T]and|9] 
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fective temperature, luminosity and mass on the a ov vs. Z plane 
were nearly unchanged in comparison with the standard models, 
but opacity enhancement had a huge impact on the instability bor- 
ders. Almost all our models had the modes £ = 0, pi and I = 1, g[ 
unstable for metallicity as low as Z = 0.007. Unfortunately, with 
the modified opacities no model reproduces the empirical values of 
/, neither for vi nor for V5 . 

As we have already mentioned, the theoretical values of the 
real part of / calculated with the standard opacities (OP, OPAL) 
agreed very well with their empirical counterparts for almost all 
modes. The problem was to fit the imaginary part of /. The change 
of the opacity near the Z-bump significantly altered the real part of 
/, whereas the imaginary part was almost unaffected. 

In the second case, opacities near the Deep Opacity Bump 
(DOB), occurring in the temperature range of log T = 6.2 - 6.5, 
were increased by 20%. This modification did not change our mod- 
els significantly. Also the instability regions for the modes V[ and 
V5 were nearly unaffected. This is caused by the fact that this 
bump is located deep inside the star and contains a small amount 
of mass. Therefore, this opacity modification did not substantially 
change the star structure. The only noticeable effect was the small 
reduction of the overshooting parameter of models fitting the /- 
parameters for both Vi and V5 modes. However, there is still no 
seismic model fitting the /-parameter for the modes V\ and v 5 si- 
multaneously. The effect of the opacity enhancement on the fre- 
quencies a nd their instabilities w ill be discussed in details in our 
next paper (Zdravkov et al.ll2013l) . 



5 CONCLUSIONS 

The aim of this paper was to give a more detailed interpretation of 
the oscillation spectrum of the hybrid pulsator 7 Pegasi. We began 
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Figure 10. The real (left panel) and imaginary (right panel) part of the nonadiabatic parameter / as a function of pulsational frequency for MODEL2 of y Peg 
marked as asterisk in Fig. 1 and Fig. 9. The empirical values of / were put for the eight well identified frequencies; circles correspond to t = 0, triangles to 
t = 1, squares to t = 2 and stars to I = 3. 



with identification of the mode degree, t, for the 14 frequencies: 
8 of the p Cep type and 6 of the SPB type. Based on the two ap- 
proaches, we were able to determine unambiguously four frequen- 
cies. For the other five frequencies two possible values of t were 
obtained. In the case of the remaining five frequencies only some 
constraints were derived. 
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Then, we tried to construct seismic models which fit the two 
low order p/g-mode frequencies and their corresponding values of 
the nonadiabatic complex parameter /. We chose the frequencies 
v\ and vj, which were identified as the radial fundamental mode 
and dipole gi mode, respectively. The problem we encountered 
was that there was no seismic model reproducing the /-parameter 
for these two frequencies simultaneously. This inconsistence can 
be caused either by the underestimated errors or/and indicate that 
some additional effects should be included in pulsation model- 
ing. One of the reasons could b e inadequacies in the opacity data. 
IZdravkov & Pamvatnvkhl (120090 have suggested increasing opacity 
by 20-50% around the Z-bump and DOB (Deep Opacity Bump) to 
explain the observed frequency range of y Peg. However, our stud- 
ies showed that these artificially increased opacities spoiled even 
more the agreement between the empirical and theoretical values 
of the /-parameter. With the modified opacities there is no model 
fitting the empirical values of /, neither for Vi nor V5. 



Although we did not fully succeeded in constructing complex 
seismic models of y Peg, we have showed directions and problems 
that need to be solved. There are the two main messages from this 
paper. The first one is a need for more accurate data on multi-colour 
time-series photometry and radial velocity data to better identify 
the observed modes and determine the empirical values of the /- 
parameter. The second one is that if opacities are to ,, blame" for 
these disagreements and problems with mode instability, the im- 
provement in computations of this microphysics data should be 
done in a more sophisticated way. 
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